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Ground state ionization potential of the He^ atom is eval- 
uated to be 5 945 204 221 (42) MHz. Along with lower order 
contributions, this result includes all eflects of the relative 
orders a*, a^me/ma and In^ a. 
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In contrast to the theoretical description of electro- 
magnetically bound two-body systems like hydrogen, 
positronium or muonium, where considerable progress is 
achieved (for the recent reviews see, e.g., |^]), high preci- 
sion calculations in more complex atoms are elaborated 
to a lesser degree. The central problem with an exten- 
sion of the methods developed for the two-body problem 
to few-electron atoms is that those methods are usually 
strongly rely on the solution of the Schrodinger equa- 
tion for a single particle in the Coulomb field. Having 
a simple analytic form, this solution is a perfect refer- 
ence point for the calculation of various observables as 
power series in the fine structure constant a using the 
quantum-mechanical perturbation theory. In higher or- 
ders of this perturbation theory, where the nonrelativistic 
approximation usually breaks down, the explicit form of 
the nonrelativistic solution facilitates the extraction of 
the ultraviolet divergences. These divergences are can- 
celed by matching to their finite counterparts calculated 
in the fully relativistic framework of the quantum elec- 
trodynamics. 

Although the Schrodinger equation for, e.g., three par- 
ticles bound by the Coulomb potentials can be solved 
numerically with very high accuracy the lack of an 
analytic solution makes the problem of the divergences 
cancellation more involved as compared to the two-body 
case. This problem was recently analyzed in Ref. |^ using 
singlet states of the helium atom as an example. Employ- 
ing the nonrelativistic quantum electrodynamics Q| reg- 
ularized dimensionally, it is demonstrated in |^ how all 
the divergences arising in the quantum-mechanical per- 
turbation theory can be extracted and canceled at the 
operator level, without recourse to an explicit form of 
the helium wave function. For the first time the 0{a^) 
correction to singlet 5* levels of the He^ atom is repre- 
sented as a sum of apparently finite average values of the 



regularization- independent operators. 

In this Letter we present the most precise evaluation 
of the helium ground state energy. Expressed in terms of 
the ionization potential (the difference between ground 
state energies of the singly charged ion and of the atom), 
our result reads: 



vth{l^S) = 5 945 204 221 (42) MHz. 



(1) 



Along with the nonrelativistic energy, this result in- 
cludes aU 0(a2), 0{a^), 0{a^) and 0{a^\n^a) rela- 
tivistic and radiative corrections. We take into account 
the finite nucleus-to-electron mass ratio M = ma/rrie = 
7 294.299 508(16) |] exactly in the nonrelativistic and 
0{a^) contributions, include the first (~ l/Af) recoil 
correction into the 0{a'^) contribution and neglect the 
nucleus recoil in higher orders. The effect of a finite nu- 
cleus charge radius i?N = 1.673(1) fm ||^ is included into 
the helium ground state energy as 
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Here Z = 2 is the nucleus charge (in units of the proton 
one), while ri and r2 denote the positions of the electrons 
with respect to the nucleus. The angle brackets in (||) and 
below denote the average value over the nonrelativistic 
ground state. In (|^) we take one half of the 0{a^ In^ a) 
correction as an estimate of the uncertainty due to higher 
orders. Our result agrees with the previous theoretical 
estimate 



DM, 
^th 1 



l^S) = 5 945 204 226 (91) MHz, 



(3) 



obtained in and including the 0{a^) and 0{a'^/M) 
effects only partially^. 

In the remaining part of this Letter we briefly describe 
details of our calculation. Ground state energy of the 
helium atom is calculated as power series in the fine 
structure constant a. Leading ('^-^ 1) contribution^, the 
Schrodinger energy E, and the corresponding wave func- 
tion ^ are found as a solution of the variational problem 



*e-mail: korobov@thsunl.jinr.su 
^e-mail: yelkhovsky@inp.nsk.su 



^The uncertainty in (^ exceeds that in (^) because only part 
of the 0(0*) corrections has been included into the calculation 
that lead to Eq.(||). 

■^Unless otherwise specified, we use the atomic units e = h = 
me = 1 and c = 1/q throughout this paper. In particular, 
the unit of energy is 2Ry = meC^Q^. 
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for the helium atom Hamiltonian taken in the nonrela- 
tivistic approximation, 
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(5) 



Here ri_2 = \ fi.2\ and r — |ri — r2|; Pi 2 are the momenta 
of the electrons and P = —pi — P2 is the momentum of 
the nucleus. 

To construct the variational wave function we use the 
simplest form of the basis, 

i;n = eM-kin-k^r2-k^r), n=l,...,N. (6) 

The complex exponents are chosen in a quasi-random 
manner from a rectangular area on the complex plane, 
for example. 



Re fc" = 



n{n + 1) 



/Pa 



Kr), (7) 



where \_x\ denotes a fractional part of x, pa is some 
prime number, while {K™™,K^^^^ is a variational in- 
terval. Imaginary parts of the parameters are generated 
in a similar way. We use both real and imaginary parts 
of ip„ to form a set of real basis functions. In particular, 
the ground state wave function ip is a linear combina- 
tion of 2N basis functions Re tpm Im ipm n = Ij ■ • ■ j 
symmetric over the interchange of the electrons positions, 
ri <-> r2. 

Variational expansion in the basis was shown in |^] 
to be very effective. It yields the best available nonrela- 
tivistic energies for many atomic and molecular systems 
and in particular for the ground state of the helium atom. 
Simplicity of the basis (|6|) allows us to evaluate analyt- 
ically matrix elements of all the operators that appear 
in the calculation. By a proper differentiation and/or 
integration of the basic integral. 



exp(-fciri - k2r2 - k^r) 



rir2r 
16^2 



(fcl -f fc2)(fc2 +A:3)(fc3 + fci) 



(8) 



with respect to fci, ^2 and k^ we express the matrix ele- 
ment of any operator involved in our calculation in terms 
of rational functions of fc's, their logarithms and diloga- 
rithms. 

For the zeroth order approximation a wave function 
built within a set of 27V = 1200 basis functions has been 
used that yields the nonrclativistic energy 



E = -2.903 304 557 727 940 23(1). 



(9) 



Here and below we cite the uncertainty of the numer- 
ical results due to finiteness of the basis set. The un- 
certainties due to incomplete knowledge of the physical 



constants are included into the final result for the ion- 
ization potential (see Table). High accuracy of (^ is not 
redundant since the calculation of rather singular matrix 
elements of higher order corrections requires very accu- 
rate variational wave function. 

First relativistic correction to the nonrclativistic value 



(y) is the average of the Breit Hamiltonian (see, e.g., 
over ip: 
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&{ri) + 5{r2) 



+ - ]: \P1-P2 + (Pi")-("P2) 

Z \ r r 

^Wl {P^k^ + (Pini)^(niP) + (1^2) 



= -1.952 050 77(1)^2 



(10) 



Here n = r/r and ni 2 = ''"i,2/'"i.2- To simplify the 
presentation, we explicitly take into consideration that 
the spin of the nucleus and the total spin of electrons are 
both equal to zero. In particular, we replace the product 
of the electron spin operators S1S2 by its eigenvalue in 
the singlet state, —3/4. 

Order and o? jM corrections to the energy can be 
represented as follows (see [|j and references therein): 
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In ri -|- 7 



miPi + (1^2; 



57.270 34(2)^3 



(11) 



Here 7 = 0.5772 ... is the Eulcr constant and ji is the 
helium Bethe logarithm |9j defined as 



/3 



_ ((pi +P2)(g - g)ln[2(g- g)](pi +p,)) 

((Pl+P2)(^-^)(Pl+P2)> 

= 4.370 039(2). 



(12) 



The cited value of /3 has been calculated for the finite 
mass of the nucleus. Details of calculations in the limit of 
no recoil (Af — s- 00) can be found in ]To| . For convenience 
of comparison with earlier results it is worth to write 
explicitly the relation to the Q-term introduced by Araki 
and Sucher 



Q = lim 



e(r - p) 

1 / In r + 7 
2^ 



(Inp + 7)(5(r) 
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The next, 0{a'^) correction to the energy is 
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I (2CnC + c2 - ^ - [J(ri) + 5{r2)] 



+ 7rfcoN('5(ri) + <5(r2)) + T:K^{5{r)) 

3C1C2CN cCnC Cnc[piP2 + "■("Pi)P2] 
4 2 4 

P^CnpI , PicVi + P2C^P2 , (PiXP2)c(PiXP2) 

8 + 4 

plc(np^f + (piw)^cp| - 3(piw)^c(np2)^ 
8 
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4 
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139.60(1)0" 



(14) 



km = - In 



3285 335 29 In 2 15C(3) 



216 547r2 
10.37657. 
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(18) 



In (|lj), all momentum operators standing to the right 
(left) of position-dependent operators are assumed to act 
on the right (left) wave function. 

Average values for a part of the operators entering into 
Eq.(p^ can be found in the literature. For the average 
values of the new operators we have obtained the follow- 
ing results: 

(VpGVp) = -0.392; (VsGVs) = -18.48; 



Cnc[PiP2 + n{np^)p2] 



= 0.811; 



pICnpI \ 
4 / 



_36.983; ( El^PllP^^ ) = 1.U2: 



= 1.078; 



(Pl X P2)c(Pl X P2) 



plc{np2Y+ {pinfcpl - 3{pinyc{np2) 



2{np2){EiP2) + {nEi) [{np2f - pI] 



-0.03; 



= -4.749; 



Here we use the following notations: C = Cn + c, 
Cn = Ci + C2, c = 1/r, Ci,2 - -^/ri,2, e = n/r^ 
and Ei,2 — — ^'T'1,2/?'i,2- The terms {VpGVp) and 
(ysGVs) in (|lj), where G is the reduced Green function 
of the Schrodinger equation, {H — E)G(rx,r2\r'i,r'2) = 
tp{i^i,'r'2)fp{'r''i,r2) — S{ri — r[)d{r2 — r'2), represent the 
effects of virtual transitions into triplet P and singlet S 
excited states, respectively (see Q for details). Pertur- 
bations which induce those transitions are 



Vp = ^1-^2 fZh _Zl2^ rxP 
4 \ rf 



where Zi_2 = T1.2 x 2; ^md 

Cn + 2c {pI + c,pI + c} Cnc 3c2 



(15) 
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Pi (Cn - c) Pi -I- (1 ^ 2) P1CP2 + (Pin)c(np2) 



(16) 



4 2 
The contact terms enter into Eq. ( |l4| ) with the coefficients 

427^2 lOZ 9ZC(3) 2179Z 
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27 47r2 
ln2 ^ 16.3557, 
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(17) 



3E1E2 - {nEi){nE2) - 2{Ei - E2)e\ 
r z > = 1.434; 



3 P2 _ 3(nP)2 
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/97rJ(r)p2\ 
-0.353; ( ) = 1.325; 



7r^(r)CN 



-2.508; 0.409; 



^ 3^J(ri)p2 



(1^2)) = 18.421; 



37r,5(ri)C2 

2 ^ 



2) 



-25.066; 



-0.958; 



2r3 V^"^"2 



(miPi Z) + (1 ^ 2) ) = 4.706. 
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Finally, the enhanced by In^ a (and hence presumably 
the leading) part of the 0{a^) correction is fl^ : 

6^'^^E = -4Z^a'^\n^{Za){6{ri) + S{r2)) « 2070a^ (19) 

Numerical results for all the contributions to the he- 
lium ionization potential are collected in the Table. Ap- 
propriate expression for the ground state energy of the 
helium ion is 
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Here /3h = 2.9841285557655 ... is the Bethe logarithm 
for the hydrogen ground state and = 3.162(2) • 10^^ 
is the nucleus charge radius i?N expressed in the atomic 
units. 

A comparison of our result (Q) and the most recent 
experimental values. 



,1S-2P 



(1^5) = 5 945 204 238 (45) MHz, 



and 



yl^-'^'^{l^S) = 5 945 204 356 (48) MHz, 



(21) 



(22) 



extracted from the measurements of 1^5* — 2^P ||l^ and 
1^5* — 2^5* [|lj] intervals, respectively, shows that the the- 
oretical value (|l|) agrees well with the former ( ^l]) and 
is within 2a from the latter ( p2[ ) if the theoretical and 
experimental uncertainties are added linearly. Efforts in 
both theoretical and experimental directions are desir- 
able in order to further clarify the situation. 

TABLE I. Contributions to the total ionization potential 
of the helium ground state. Uncertainty in the nonrelativistic 
value is due to uncertainty in the nucleus mass. 

^^'th(l'S'), MHz 

5 945 262 288.62(4) 
-16 800.338(4) 
-40 486.375(50) 
-834.9(2) 
84 (42) 
-29.55(4) 
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